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TECHNICAL  NOTE  3508 


LAMINAR  FREE  CONVECTION  ON  A  VERTICAL  PLATE  WITH  PRESCRIBED  NONUNIFORM 
WALL  BEAT  FLUX  OR  PRESCRIBED  NONUNIFORM  WALL  TEMPERATURE 

By  E.  M.  Sparrow 


SUMMARY 

An  analysis  is  made  for  laminar  free  convection  on  a  vertical  plate 
with  nonuniform  thermal  conditions  at  the  surface.  Prescribed  variations 
are  considered  for  the  wall  heat  flux  and  for  the  wall  temperature. 

For  the  situation  where  the  wall -heat -flux  variation  is  prescribed, 
graphs  are  presented  from  which  the  resulting  wall -temperature  variation 
may  be  obtained.  Local  heat-transfer  coefficients  may  be  readily  deter¬ 
mined  using  the  information  given  on  the  graphs.  Results  for  the  impor¬ 
tant  special  case  of  uniform  wall  heat  flux  are  also  given. 

For  the  situation  where  the  wall -temperature  variation  is  prescribed, 
graphs  are  presented  from  which  the  over-all  rate  of  heat  transfer  from 
any  length  of  the  plate  may  be  obtained.  Another  set  of  graphs  is  pre¬ 
sented  for  obtaining  local,  heat-transfer  coefficients . 

All  the  aforementioned  results  are  given  for  fluids  having  Prandtl 
numbers  in  the  range  0.01  to  1000. 

The  flow  is  taken  to  be  of  the  boundary-layer  type,  and  the  problem 
is  formulated  by  the  Karman -Pohlhaus en  method.  The  solution  of  the  result¬ 
ing  equations  is  achieved  by  series  expansion.  The  first  term  of  the 
series  corresponds  to  the  result  for  uniform  thermal  conditions  on  the 
wall.  The  succeeding  terms  give  the  influence  of  the  nonuniform  thermal 
conditions.  The  first  five  terms  of  the  series  have  been  calculated. 


INTRODUUTION 

T.»minftT  free  convection  on  a  vertical  plate  has  been  a  subject  of 
study  since  1881.  Most  of  the  analytical  work  has  been  done  for  the  sit¬ 
uation  where  the  wall  temperature  is  uniform  over  the  entire  surface. 

An  exact  solution  of  the  boundary-layer  differential  equations  for  free 
convection  on  a  vertical  flat  plate  with  uniform  wall  tenperature  is 
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given  In  reference  1  for  several  Prandtl  numbers  In  the  range  0.01  to 
1000.  Reference  1  also  lists  other  work  on  the  uniform- wal  1. -temperature 
problem.  The  Karmen -Pohlhaus en  method  is  applied  to  the  uniform-wall- 
temperature  case  in  references  2  (pp.  671-673)  and  3  (pp.  158-163).  The 
results  reported  in  these  references  agree  well  with  those  from  the  ex¬ 
act  solution  of  reference  1. 

Considerably  less  analytical  work  has  been  done  for  the  situation 
where  the  heat  transfer  is  uniform  over  the  surface.  An  exact  solution 
of  the  boundary-layer  differential  equations  for  the  uniform-heat -flux 
case  is  given  in  reference  4,  which  also  contains  results  calculated  by 
the  Karman -Pohlhaus  en  method.  The  results  of  the  exact  solution  and 
those  from  the  KR-rrnan -Pnhl hanssn  method  are  in  good  agreement. 

Accounts  of  experimental  investigations  of  free  convection  on  ver¬ 
tical  surfaces  are  given  in  references  1  to  6. 

In  a  large  number  of  technical  applications  the  thermal  conditions 
on  the  surface  are  nonuniform.  '  These  nonuniformities  in  thermal  condi¬ 
tions  may  be  grouped  into  two  categories: 

(1)  The  heat  flux  may  be  prescribed  to  vary  over  the  surface.  It 
is  then  of  interest  to  calculate  the  resulting  variation  of  the  surface 
temperature . 

(2)  The  variation  c£  the  temperature  on  the  surface  may  be  pre¬ 
scribed.  It  is  then  of  interest  to  calculate  either  the  local  rate  of 
heat  transfer  at  various  locations  on  the  surface,  or  the  over-all  rate 
of  heat  transfer  from  the  surface,  or  both. 

This  report  presents  a  first  attempt  at  solution  of  the  free- 
convection  problem  on  a  flat  plate  for  these  two  categories  of  non- 
uniform  thermal  conditions  at  the  surface. 

This  analysis  was  made  at  the  MCA  Lewis  laboratory. 


GENERAL  CONSIDERATIONS 

The  physical  problem  and  the  coordinate  system  are  indicated  in  the 
following  sketches,  which  show  a  vertical  surface  that  may  represent  a 
flat  plate  or  a  vertical  cylinder  of  large  diameter: 
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Sketch  (a).  -  Heat  transfer  from  Sketch  (b).  -  Heat  transfer  from 
wall  to  fluid  and  t  >  t  fluid  to  wall  and  t  <  t 

everywhere .  everywhere  w  a 

Sketch  (a)  depicts  a  situation  in  which  the  heat  transfer  is  from 
the  wall  to  the  fluid  at  all  points  on  the  surface  and  the  local  wall 
temperature  tw  everywhere  exceeds  the  ambient  temperature  t&.  On  the 
wall.,  either  the  heat  flux  may  be  designated  to  vary  with  x,  or  the 
wall  temperature  may  be  prescribed  to  vary  with  x.  (No  variations  are 
considered  in  the  direction  of  the  z-coordinate,  normal  to  the  pagej 
the  problem  is  thus  taken  to  be  two  dimensional.)  The  fluid  in  the 
neighborhood  of  the  wall  has  a  higher  temperature  and  a  lower  density1 
than  the  fluid  far  from  the  wall.  Thus,  because  of  buoyancy,  there  will 
be  established  an  upward  flow  of  fluid  in  the  neighborhood  of  the  wall. 
The  region  of  space  in  which  the  upward  flow  primarily  occurs  is  called 
the  velocity  boundary  layer.  A  thermal  boundary  layer  is  defined  as  that 
region  of  space  where  the  temperature  t  deviates  markedly  from  the  am¬ 
bient  temperature  t&.  In  general,  the  velocity  and  thermal  boundary 

layers  have  different  thicknesses,  the  relative  magnitudes  depending  upon 
the  fluid  properties .  Both  boundary  layers  are  assumed  to  have  zero 
thickness  at  the  leading  edge  (x  =  0) .  The  velocity  boundary  layer  is 
shorn  schematically  in  the  sketch. 

Sketch  (b)  shows  a  situation  in  which  the  heat  transfer  is  from  the 
fluid  to  the  wall  at  all  points  of  the  surface  and  the  local  wall  temper¬ 
ature  tw  is  everywhere  less  than  t&.  Again,  on  the  wall,  either  the 
heat  flux  or  the  temperature  may  be  prescribed  to  vary  with  x.  Here, 
the  flow  of  fluid  in  the  boundary  layer  is  downward  as  shown. 


1This  refers  to  fluids  showing  the  usual  trend  of  density  decreasing 
with  increasing  temperature. 
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If  the  coordinate  systems  are  taken  as  shown  in  the  sketches,  the 
method  of  analysis  and  the  results  for  the  heat-transfer  parameters  are 
the  same  for  these  two  situations;  and  there  will  be  no  need  to  treat 
them  separately.  So,  the  analysis  will  be  carried  out  for  the  case  of 
heat  transfer  from  the  wall  to  the  fluid  (t  >  t&) ,  but  it  is  to  be  remem¬ 
bered  that  the  results  apply  to  both  situations  depicted  in  the  sketches. 

A  definite  class  of  variations  of  the  wall  heat  flux  wi.ll  be  pre¬ 
scribed'.  Suppose  that  the  region  of  interest  on  the  vertical  plate  lies 
between  x  =  0  and  another  location  x  =  x^  (x-^  must  lie  in  the  region 

of  laminar  flow  over  the  plate).  The  wall  heat  fluxes  to  be  considered 
here  have  a  finite  nonzero  value  at  x  =  0  and  either  increase  steadily 
from  x  =  0  to  x  =  x^  or  else  decrease  steadily  from  x  =  0  to 
x  =  xL. 

The  form  of  the  wall -temperature  variations  considered  here  is  sim¬ 
ilar  to  that  outlined  for  the  wall -heat -flux  variations.  The  wall  tem¬ 
perature  relative  to  ambient  has  a  finite  nonzero  value  at  x  =  0  and 
either  rises  steadily  from,  x  =  0  to  x  =  x-^  or  else  decreases  stead¬ 
ily  from  x  =  0  to  x  =  x^. 

Although  the  analysis  made  here  is  for  free  convection  in  a  gravity 
field,  it  may  easily  be  generalized  to  include  other  force  fields.  It 
3s  only  necessary  to  replane  the  gravitational  force  per  unit  mass  g 
by  the  body  force  per  unit  mass  of  the  other  force  field  under 
cons iderat ion . 


BASIC  EQUATIONS 


The  equations  expressing  conservation  of  mass,  momentum,  and  energy 
for  steady  laminar  flow  in  a  boundary  layer  on  a  vertical  flat  plate  are 
as  follows: 


du  dv  _  n 


c)u  du  n  ( , 

u  TI  +  v  x*  =  SPCb 

ox  ay 


ta)  +  v 


&t  St 
u  —  +  v  —  =  a 


dx 


dy 


^t 

By2 


(i) 

afu 

(2) 

Sy2 

(3) 

(AH  symbols  are  defined  in  appendix  A.)  In  accordance  with  the  usual 
practice  in  free  convection,  the  density  is  considered  a  variable  in 
formulating  the  buoyancy  term  gp(t  -  ta) .  Aside  from  this,  the  fluid 

properties  are  taken  constant.  Viscous  dissipation  and  work  against 
the  gravity  field  are  neglected. 
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Following  reference  3,  It  is  assumed  that  a  common  boundary-layer 
thickness  6  can  be  used  for  both  velocity  and  thermal  boundary  layers. 
This  assumption  has  its  justification  in  the  fact  that  the  results  of 
calculations  performed  with  it  are  in  good  agreement  with  those  from 
exact  solutions  of  the  boundary-layer  differential  equations  for  the 
cases  of  uniform  wall  temperature  and  uniform  heat  flux. 


Then,  equations  (l)  to  (3)  are  integrated  across  the  boundary 
layer  to  give 


_d_ 

dx 


u(t 


(3a) 


The  integrated  form  of  equation  (l)  has  been  absorbed  into  equations  (2a) 
and  (3a) . 

These  equations  have  a  definite  physical  meaning.  They  are,  in  fact, 
expressions  of  the  conservation  laws  for  the  element  of  boundary  layer 
shown  in  the  following  sketch: 


Sketch  (c) 

Solutions  for  the  velocity  and  temperature  in  the  boundary  layer 
will  be  obtained  which  satisfy  these  conservation  equations  and  the 
boundary  conditions.  These  solutions  will  in  turn  be  vised  to  calculate 
the  important  heat -transfer  parameters. 
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ANALYSIS 


Prescribed  Nonuniform.  Nall  Heat  Flux 


The  analysis  proceeds  ■with  the  use  of  the  Karman-Pohl  hausen  method, 
according  to  which  the  velocity  and  temperature  distributions  in  the 
boundary  layer  are  written  as  polynomials  in  y  whose  coefficients  are 
functions  of  x.  The  coefficients  are  found  from  the  boundary  conditions 
of  the  problem  and  by  using  the  Integrated  momentum  and  energy  equations 
(eqs.  (2a)  and  (3a)). 


The  following  polynomials  are  chosen: 


u 


=“IM)2 


(4) 

(5) 


Equation  (4)  satisfies  the  conditions: 
<dt 


¥ 


=  -  §  when  y  =  0;  t  =  t 


and  ^  =  0  (smooth-fit  condition)  when  y  =  5.  The  wall  heat  flux  q 
is  to  be  regarded  as  a  specified  function  of  x.  The  conditions  satis¬ 
fied  by  equation  (5)  are:  u  =  0  when  y  =  0;  u  =  0  and  =  0 

(smooth-fit  condition)  when  y  =  5.  The  functions  ca(x)  and  6(x)  still 
remain  to  he  determined. 

The  polynomials  representing  the  velocity  and  temperature  distribu¬ 
tions  are  introduced  into  equations  (2a)  and  (3a),  and  after  the  inte¬ 
gration  is  carried  out,  there  results  a  pair  of  first  order,  ordinary 
differential  equations  for  o>  and  8.  In  dimensionless  form,  these 
equations  are 


105  dX 


(22A)  =  %- 


ffi 

A 


(6) 


30  dX  (2A 


2  2_  _q_ 

oJ 


%)  % 


(7) 


X,  2,  and  A  are  the  dimension] ess  counterparts  of  x,  cd,  and  8  and 
axe  defined  in  the  Symbol  list  (appendix  A) .  The  symbol  qo,  which 
represents  the  heat  flux  at  the  leading  edge  (x  «=  0) ,  is  used  as  a 
reference  heat  flux  in  the  rest  of  the  analysis. 
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Next,  the  variation  of  q/c^  is  to  be  prescribed.  Suppose  that 
the  region  of  interest  on  the  vertical  plate  lies  between  x  <=  0  and 
another  location  x  =  x^.  The  surface  heat  fluxes  to  be  considered  here 
have  finite  nonzero  values  at  x  =  0  and  either  increase  steadily  from 
x  =  0  to  x  =  x-^  or  else  decrease  steadily  from  x  =  0  to  x  =  x^. 
Explicitly,  the  class  of  surface-heat-flux  variations  to  be  considered 
is  written  as 


t =  1  *  £(tf  (8) 

The  exponent  r,  which  gives  the  shape  of  the  variation,  is  required  to 
be  a  positive  number,  integral  or  nonintegral.  From  equation  (8),  it 
may  be  seen  that  s  represents  the  maximum  (percentage)  deviation  of 
q  from  q^  in  the  region  x  =  0  to  x  =?  xL. 

For  the  q/c^  given  by  equation  (8),  the  pair  of  differential 
equations  (6)  and  (7)  can  be  solved  by  expanding  2  and  A  in  Maclaurin 
series  in  terms  of  6  as  follows: 

2(x,e,Pr,r)  =  ^  e^folh^r)  =  Qq^Pt^) 
n=0 

CO 

A(X,e,Pr,r)  =  ^  e^Cx^r)  =  ^(X,Fr,r) 
n=0 

The  expressions  for  2,  A  and  q/qQ  are  introduced  into  equations 

(6)  and  (7).  Terms  are  grouped  according  to  the  power  of  e  that  mul¬ 
tiplies  them,  that  is, 

]♦--<>  C-> 

-s]  +  el[  ]+ “°  Cl2) 

In  order  that  equations  (ll)  and  (l2)  be  satisfied  for  any  value  of  6, 
each  of  the  brackets  must’ be  identically  zero.  Equating  to  zero  the 
brackets  multiplying  e®  yields  a  pair  of  simultaneous  equations  for 
2q  and  A^.  In  a  similar  fashion,  the  brackets  multiplying  e11  give 
a  pair  of  equations  for  2n  and  A^.  It  may  be  noted  that  the  equations 
for  2n  and  Aq  will  include  all  the  20  through  2q_1  and  the  Aq 
through  . 


+  c21(X,Pr,r)  +  (9) 

+  cA1(X,Pr,r)  +  •  •  •  (lO) 
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The  solution  for  2q  and  Aq,  which  coincides  with  the  results 
for  uniform  wall  heat  flux  (e  =  0),  is  treated  more  fully  in  appendix 
B.  The  succeeding  2n  and  AQ  evidently  give  the  effects  of  the  non- 
uniformity  of  the  heat  flux. 

The  first  five  functions  in  the  Maclaurin  series  (2q  through  2^ 
and  Aq  through  A^_)  have  "been  computed.  The  results  are  listed  in 
concise  form  as  follows: 


\  -  bbn  xn(Pr’r)^; 

\1XC  / 

9^ 

(13) 

S2n  -  (±l)“  rn( Fr,r/j 

L^x3/5 

(14) 

The  factors  Xn  and  were  calculated  from  linear  algebraic  equations. 

The  results  thus  found  are  applied  in  a  later  section  in  which  the 
heat -transfer  parameters  are  obtained. 


Prescribed  Nonuniform  Wall  Temperature 

The  analysis  for  the  pres cribed-nonuniform-wall- temper ature  case 
proceeds  in  a  similar  fashion  to  that  for  the  prescribed  heat-flux  case. 
The  temperature  and  velocity  profiles  in  the  boundary  layer  are  approxi¬ 
mated  by  the  following  polynomials : 

-  I)  s  K1  - 1)  C15> 

u  “  ®  I  (x  - 1)  be) 

Equation  (l5)  satisfies  the  conditions:  t  =  t  when  y  =  0;  t  =  ta 
and  =  0  when  y  =  8.  The  wall  temperature  t^  is  to  be  regarded 

as  a  specified  function  of  x.  Since  the  ambient  temperature  ta  is 
taken  to  be  constant,  then  9  is  regarded  as  a  prescribed  function  of 
x.  The  conditions  satisfied  by  equation  (l6)  are:  u  =  0  when  y  =  0; 

u  =  0  and  =  0  when  y  =  5. 

oy 

The  variation  of  the  wall  tengperature  t^.  is  prescribed  in  a  re¬ 
gion  of  interest  between  x  =  0  and  some  other  location  x  =  x^.  The 

wall  temperature  relative  to  ambient  considered  here  has  a  finite  non¬ 
zero  value  at  x  =  0  and  will  either  increase  steadily  from  x  =  0  to 
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x  =  x^  or  else  decrease  steadily  from  x  =  0  to  x  =  x^.  The  explicit 
form  of  the  temperature  variation  considered  here  is 


Again,  the  exponent  r  is  a  positive  number,  integral  or  nonintegral. 

As  may  be  seen  from  equation  (l7a),  e  is  the  maximum  (percentage)  de¬ 
viation  of  0  from  0Q  in  the  region  0  <  x  <  x^. 

Equations  (l5)  and  (l6)  are  introduced  Into  the  integrated  momentum 
and  energy  equations  (eqs.  (2a)  and  (3a)).  The  resulting  simultaneous 
equations  for  a)  and  5  are  solved  subject  to  the  prescribed  wall- 
temperature  variation  (eq.  (l7))  using  Maclaurin  series.  The  first  five 
functions  in  the  Maclaurin  series  have  been  calculated.  The  results 
have  a  form  similar  to  those  given  in  equations  (l3)  and  (l4).  The 
solution  for  the  uniform-waU-temperature  case  coincides  with  the  lead¬ 
ing  term  of  the  Maclaurin  series  expansion. 


HEAT -TRANSFER  RESULTS 
Prescribed  Nonuniform  Wall  Heat  Flux 


For  the  case  where  the  wall  heat  flux  is  prescribed,  it  is  of  in¬ 
terest  to  determine  the  resulting  wall -temperature  variation  and  local 
heat -transfer  coefficients . 


Resulting  wall -temperature  variation.  -  The  wall  temperature  is  ob¬ 
tained  by  setting  y  equal  to  0  in  equation  (i).  Thus, 


t  -  t 

\  “  2k 


(18) 


or 


w 


2k 


2k 


(18a) 


For  a  uniform  wall  heat  flux  q  =  qQ,  it  has  already  been  noted  that 
A  =  dQ}  so, 


- 

'  w 


ta)q 


2k 


(19) 
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■where  (t  -  t  )  is  the  wall -temperature  variation  corresponding  to  a 

W  a  Hq 

■uniform  heat  flux  q  =  q^.  In  equation  (l8a) ,  A/A^  is  evaluated  from 
equations  (8),  (10),  and  (l3),  and  qQ8o/2k  is  evaluated  from  equation 
(l9),  giving  the  result 


(20) 


2 

This  equation  gives  the  ratio  of  the  wall  temperature  at  some  loca¬ 
tion  x  (x  <  xL)  having  a  specified  ratio  i/Iq  to  the  wall  temperature 

at  the  same  location  on  a  plate  having  a  uniform  heat  flux  q  =  qQ. 

Since  the  X's  depend  upon  the  Prandtl  number  and  r,  the  temperature 
ratio  given  by  equation  \20)  depends  on  q/qQ,  Pr  and  r. 


The  ratio 


(t. 


w 


*a) 


is  plotted  in  figures  l(a)  to  (e) .  Each  of 


the  plots  applies  for  a  specific  value  of  r.  The  values  of  r  =  0,  l/2, 
1,  2,  and  3  have  been  used  for  the  five  plots.  Results  for  other  values 
of  r  may  be  obtained  by  replotting''’  the  information  given  in  figures 
l(a)  to  (e)  using  r  as  the  abscissa  variable.  On  each  of  the  plots, 
the  temperature  ratio  is  plotted  against  q/qQ  for  the  range 


0.5  <q/qQ<  1.5  with  Prandtl  numbers  between  0.01  and  1000  appearing 
as  parameters  on  the  curves.  In  cases  where  curves  for  different  Prandtl 
numbers  fall  so  close  together  as  to  make  it  impossible  to  plot  them 
separately,  one  curve  was  used  for  the  several  Prandtl  numbers. 


Once  the  ratio  y  - — r-y —  has  been  determined  from  one  of  the 
l>tw  “  ta'q0 

figures  (or  by  replotting  the  data  given  therein),  the  wall  temperature 
(tw  -  ta)  can  be  found  if  an  expression  for  (tw  -  ta)^  is  given.  The 


expression  for  (tw  -  ta)^ 

cane  of  uniform  heat  flux 
Kerman -Pnhl hausen  method. 


,  the  wall  temperature  corresponding  to  the 

q  =  qQ,  is  derived  in  appendix  B  by  the 
and  appears  on  figure  1. 


“^Relative  to  ambient"  is  to  be  understood  every  time  the  wall  tem¬ 
perature  is  mentioned. 

^r  =  0  is  a  limiting  case,  which  is  included  here  to  facilitate 
the  replotting  when  results  for  r  near  zero  are  required. 
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Local  heat-transfer  coefficient.  -  The  definition  of  the  local 


OJ 

CM 

o 


heat-transfer  coefficient  is 


hx  t  -  t 

8. 

For  the  case  of  a  uniform  heat  flux  q  =  q  ,  equation  (21)  is 


**>*<>  "  TV  -  *a) 


(21) 


(21a) 


The  ratio  of  equation  (2l)  to  equation  (21a)  is 

kx  _  _q_  ^ 

Sd  (tw  “  ta^ 


(22) 


This  equation  gives  the  ratio  of  the  local  coefficient  at  some  position 
x  having  a  specified  ratio  q/qQ  to  the  local  coefficient  at  the  same 
position  on  a  plate  having  a  uniform  wall  heat  flux  q  =  qQ-  When  q/q^, 

Pr  and  r  are  known,  the  temperature  ratio  appearing  in  equation  (22) 
is  known  (from  fig.  l)j  hence,  h./h  can  he  calculated.  The  expres- 

3  ^o 

sion  h_  n  is  given  in  appendix  B. 

‘Prescribed  Nonuniform  Wall  Temperature 

For  the  case  where  the  wall  temperature  iB  specified,  it  is  of  in¬ 
terest  to  calculate  the  over-all  heat -transfer  rate  and  local  heat- 
transfer  coefficients . 


Over-all  heat-transfer  rate.  -  The  over-all  rate  of  heat  transfer 
in  a  section  of  plate  of  width  b  from  x  =  0  to  another  location 
x(x  $  x-jO  is 

Q  =  b  q  dx  =  -  kb C23) 


= b  r 

MD 


The  derivative  (dt/§y)^_Q  is  evaluated  from  equation  (l5) .  Substitution 
into  equation  (23)  gives 


Q  =  2kb  /  ^  dx  =  2kb 

Jo  5 


Xi 


dx 


S 


(24) 
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For  the  case  of  uniform  wall  temperature  9  =  9Q ,  equation  (24)  hecomes 

The  ratio  of  equation  (24)  to  equation  (25)  is 


Equation  (26)  is  the  ratio  of  the  over-all  heat  transfer  from 
x  =  0  to  x  (x  <  x^)  on  a  plate  with  variable  surface  temperature  to 
the  heat  transferred  in  the  same  region  on  a  plate  with  a  uniform  wall 
temperature  9  -  9  .  Since  8  corresponding  to  the  prescribed  wall- 
temperature  variation  (eq.  (l7))  has  been  calculated  in  the  Analysis 
section,  it  is  possible  to  evaluate  the  integrals  appearing  in  equation 
(26). 

The  results  for  q/q^  ,  which  are  shown  in  figures  2(a)  to  (e), 

depend  upon  the  value  of  0/0Q  at  the  location  x,  upon  Prandtl  num¬ 
ber,  and  upon  r.  Each  plot  is  for  a  specific  value  of  r.  The  values 
of  0,  l/2,  1,  2,  and  3  have  been  used  for  the  five  plots.  Results  for 
other  values  of  r  may  be  obtained  by  replotting  the  information  given 
in  figures  2(a)  to  (e)  with  r  as  the  abscissa  variable.  On  each  plot, 
the  ratio  Q/Q q  is  plotted  against  0/0o  for  the  range  0.7< 0/@o<  1.3, 

with  Prandtl  numbers  between  0.01  and  1000  appearing  as  a  parameter  on 
the  curves . 

It  is  to  be  emphasized  that  the  value  of  0/0Q  x  is  to  be 

used  for  determining  values  from  the  figure  when  the  over-all  heat  trans¬ 
fer  from  x  =  0  to  x  is  required.  The  quantity  ,  which  is  the 

o 

over-all  heat-transfer  rate  from  x  =  0  to  x  on  a  plate  having  uni¬ 
form  surface  temperature  9  =  0Q,  is  calculated  in  reference  3  by  the 

Kft-rman-Pohl hausen  method.  The  expression  for  Qg  is  given  in 
figure  2.  0 

Local,  heat-transfer  coefficient.  -  From  the  defining  equation  (eq. 
(2l)),  the  local  coefficient  is  found  to  be 
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The  derivative  is  obtained  from  equation  (l5) .  For  uniform  wall  temper¬ 
ature  0  =  5  -  V  therefore 

Bq  =^0o 

where  tu.  n  is  the  local  coefficient  for  the  case  of  a  uniform  wall 
temperature  0  =  0Q.  Combining  equations  (27)  and  (28)  gives 

% 

"  6 

This  ratio  is  plotted  in  figures  3(a)  to  (e)  in  the  manner  already  out¬ 
lined  for  the  preceding  results  of  the  analysis .  The  expression  for 
q  (from  ref.  3)  is  given  on  the  figure. 


DISCUSSION 

The  following  generalizations  can  be  made  from  inspection  of  the 
results : 

(1)  At  a  location  x  with  a  fixed  q/q.  the  deviation  of 
tw  -  ta 

-7-7 - — -t —  from  1  increases  with  decreasing  r  for  a  fixed  Prandtl 

“  ta 'qo 

number.  Also,  the -deviation  from  1  is  larger  for  larger  Prandtl  numbers 
at  a  fixed  r  and  q/q^. 

(2)  In  a  length  of  plate  between  x  =  0  and  x,  having  a  fixed 
0/0q  at  x,  it  is  seen  that  the  deviation  of  Q/Q q  from  1  increases 

with  decreasing  r  at  a  fixed  Prandtl  number.  Also,  the  deviation  from 
1  is  smaller  for  the  larger  Prandtl  numbers  at  a  fixed  r  and  0/0Q. 

(3)  The  deviation  of  tu/h  from  1  increases  with  increasing  r 

X,qQ 

at  a  fixed  Prandtl  number  and  a  fixed  q/qQ  at  x.  Also,  the  deviation 
from  1  is  smaller  for  the  larger  Prandtl  numbers  at  a  fixed  r  and  q/qQ. 

(4)  The  deviation  of  tu/h^  g  from  1  increases  with  increasing  r 

’  o 

at  a  fixed  Prandtl  number  and  a  fixed  0/0Q  at  x.  Also,  the  deviation 
from  1  is  smaller  for  the  larger  Prandtl  numbers  at  a  fixed  r  and 

8/V 
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It  has  already  been  noted  that  the  temperature  ratio 


^ 


•was  plotted  in  figure  1  for  0.5  ^  q/qQ<1.5.  This  range  was  decided 

upon  by  studying  the  extent  of  the  errors  due  to  truncation  of  the 
infinite  series  in  equation  (20)  after  n  =  4.  The  range 
0.7  <  9/Qq^  1.3  used  in  figures  2  and  3  wan  decided  upon  by  studying 
the  truncation  errors  of  series  associated  with  equations  (26)  and  (29). 


CONCLUDING  REMARKS 

When  the  nonuniform,  wall  heat  flux  can  be  written  in  the  form  of 
equation  (8),  the  resulting  wall  temperatures  may  be  found  directly  from 
figure  1,  and  the  local  heat-transfer  coefficient  can  be  calculated  from 
equation  (22).  All  that  is  needed  to  use  the  graphs  is  the  value  of 
q/qQ  at  the  point  of  interest,  the  exponent  r  which  gives  the  shape 
of  the  q/qQ  variation,  and  the  Prandtl  number.  It  is  to  be  noted  that 
e  is  not  needed. 

When  the  wall  temperature  is  specified  by  a  relation  of  the  form  of 
equation  (17),  the  over-all  heat-transfer  rate  and  local  coefficients  are 
found  from  figures  2  and  3,  respectively.  The  over-all  heat -transfer 
rate  for  a  section  of  plate  from  x  =  0  to  x  (x  xL)  may  be  found  from 

the  graphs  when  the  following  are  known:  0/0Q  at  x,  the  exponent  r 
which  gives  the  shape  of  the  0/dQ  variation,  and  the  Prandtl  number. 

The  same  quantities  are  needed  to  find  the  local  heat-transfer  coefficient 
at  x. 

It  is  recognized  that  whenever  a  new  application  of  the  Karman- 
Pohlhausen  method  is  made,  it  is  desirable  to  confirm  the  results  by 
checking  with  those  of  experiment  or  of  a  less  approximate  analysis. 

The  author  is  not  acquainted  with  any  experimental  data  that  may  be 
used  to  check  the  results  derived  here  for  nonuniform  thermal  conditions 
at  the  surface.  Nor  is  there  now  available  any  other  analysis  with  which 
the  present  results  may  be  compared. 

For  the  special  cases  of  uniform  wall  temperature  and  uniform  heat 
flux,  there  are  exact  solutions  of  the  laminar-boundary-layer  equations 
aa  well  as  experimental  data.  For  these  cases,  the  heat-transfer  re¬ 
sults  derived  from  the  Karman -Pohlhaus en  method  agree  well  with  those 
from  the  exact  solutions  and  those  of  experiment. 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  May  4,  1955 
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APPENDIX  A 

SYMBOLS 


g 


Gr 


Gr 


h 

x 


n 

Pr 

Q 


r 

t 

u 

v 

X 


XL 

X 


The  following  symbols  are  used  in  this  report: 
width  of  plate,  ft 

specific  heat  at  constant  pressure,  Btu/(lb)  (<^F) 
acceleration  due  to  gravity,  ft/see2 

gPlo*4 

modified  Grashof  number  based  on  x,  - 5—,  dimensionless 

3  kv 

Grashof  number  based  on  x,  - - — ,  dimensionless 

local  heat-transfer  coefficient,  Btu/(sec)(sq  ft)(°P) 

thermal  conductivity,  Btu/(sec) (ft) (°F) 

index  for  naming  terms  of  a  series,  dimensionless 

C  Lt 

p^  v 

Prandtl  number,  ,  dimensionless 

’  k  as 

over-all  heat-transfer  rate  on  a  plate  of  width  b  be¬ 
tween  x  =  0  and  x  (x  <  xL),  Btu/sec 

local  heat-transfer  rate  per  unit  area,  Btu/(sec)(sq  ft) 

heat-transfer  rate  per  unit  area  at  x  =  0,  Btu/(sec)(sq  ft) 

exponent  defined  by  eqs .  (8)  and  (l7),  dimensionless 

static  temperature, 

velocity  component  in  x-direction,  ft/sec 

velocity  component  in  y-direction,  ft/sec 

coordinate  measuring  distance  along  plate  from  leading 
edge,  ft 

coordinate  defining  region  of  interest  x  =  0  to  x  =  x^. 


ft 

dimensionless  coordinate. 


gpqA1/4 

tv2/ 


X 
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y 

a 


Vri*  * 


o 


0 

0 


0’ 


n 

v 

o 

a> 


'  T, 


AoA'  *  *  >Aa 


V 


2 

^0>®1>  * 


•2t 


coordinate  measuring  normal  distance  from  plate ,  ft 


thermal  diffusivity,  - ,  sq  ft/sec 


PCt 


coefficient  of  thermal  expansion,  (°F) 


-1 


factors  in  solutions  for  2q,2j_, — 2n  defined  by  eq. 
(l4),  dimensionless 

dimensionless  boundary -layer  thickness  for  prescribed- 

1/4 


heat -flux  case 


■m 


coefficients  in  Maclaurin  series  expansion  of  A  in 
terms  of  £,  dimensionless 

boundary -layer  thickness,  ft 

boundary -layer  thickness  for  uniform  heat  flux  or  uniform 
vail  temperature,  ft 

number  giving  percentage  deviation  of  q  at  x  =  x^  from 
qQ,  or  number  giving  percentage  deviation  of  0  at 
x  =  x^  from  0Q,  dimensionless 

vail-  to  ambient -temperature  difference,  t  -  ta,  °F 

vail -to  ambient- temperature  difference  at  x  =  0,  °F 


•A^  defined  by  eq. 


factors  in  solutions  for  An,A  , 

(l3),  dimensionless 

absolute  viscosity,  lb/(sec)(ft) 

kinematic  viscosity,  sq  ft/sec 

density,  lb/cu  ft 

velocity  function  defined  by  eqs.  (5)  and  (l6),  ft/sec 

(gPqo^Y1  4 

dimensionless  velocity  function,  V - g — J  CD 

coefficients  in  Maclaurin  series  expansion  of  2  in 
terms  of  6,  dimensionless 
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Subscripts: 

a  ambient 


q.  on  a  plate  having  uniform,  heat  flux 

w  wall 

0Q  on  a  plate  having  uniform  wall-  to  ambient-temperature 

difference  0Q 
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APPENDIX  B 


RESULTS  FOR  UNIFORM  WALL  HEAT  FLUX 
Men  q'  =  for  all  values  of  x,  equations  (6)  and  (7)  become 


JL  A  (A)  _  fl 

105  dX  v^cro'  “6  Aq 


1  d  /q  . 2\  _  2 

30  dX  nM(A)'  ~  P? 


(Bl) 


(B2) 


The 


solutions  for  &0  and  Aq  are 


20  =  (6000)1/5  Pr_1/5^|  +  Prj  7  X3/ 


Aq  =  (360)1//5^|  +  Prj  ^  Pr“2/5  X1/5 
Equation  (B4)  can  be  rewritten  as 

-  -I1/5 

IT  =  LT  =  (360)1/5! 


X 


(I  *  Pr)' 


P2  J 


where  Gr*  is  a  modified  Grashof  number  based  on  x  and  defined  by 


(B3) 

(B4) 

(B5) 


Gr*  = 
x 


gPV< 


kv 


(B6) 


The  surface -temperature  distribution  found  by  introducing  equation  (B5) 
into  equation  (l8)  is 


(t  -  t  )  =  1.622  f  °rj-±  g£\  7 

»  W  n  ‘  n  _ 


a  9q 


yPr2  Gr 


(B7) 


The  local  heat-transfer  coefficient  is  obtained  by  the  following 
rearrangement  of  equation  (B7):  /,_ 

qp  o.62kfPr2^x>\ 7 

^qo  ~  ^qo  X  ^°’8  +  ^ 


(B8) 


i 
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Figure  1,  -  Continued.  Ratio  of  well  temperature  at  bo»o  location  x  on  plate  with  variable  heat 
flux  to  wall  temperature  at  ea«»  location  on  plate  with  unlfom  heat  flux  q  ■  Q0  •  Absoloea  Is 
value  of  q/q.,  at  point  of  interest;  r  1b  an  exponent  In  equation  (S) (relation  oy  whloh  heat- 
flux  variation  la  preeorlbed);  Prandtl  number  1b  parameter  on  ourveB. 


Figure  1.  -  Continued.  Ratio  of  wall  taxperaturo  at  some  location  x  on  plate  with  variable  heat 
flux  to  wall  temperature  at  cam©  location  on  plate  with  uniforw  heat  flux  q  ■  q0  ♦  Abaoieea  le 
value  of  q/*^  at  point  of  interest;  r  is  an  exponent  in  equation  (a)  (relation  by  whioh  heat- 
flux  variation  in  prescribed))  Prandtl  nixiber  is  parameter  on  curves. 


(a)  V  -  3. 

K*11?  °f  ™ll  temperature  at  some  location  x  do  plata  with  variable  heat 
flux  to  wall  temperature  at  sane  location  on  plate  with  uniform  heat  flux  q  -  a  .  Abscissa  Is 

iat  POl“t,KL>ntSrS^r  lB  ““  *I**:,nent  111  aquation  (B) (relation^/  which  heat- 
flux  variation  ia  prescribed)  j  Prandtl  number  la  parameter  on  curve* . 
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(a)  r  m  0. 


Figure  2.  -  Ratio  of  over-all  heat  transfer  from  x  ■»  0  to  x  on  plai 
with  variable  wall  temperature  to  over-all  heat  transfer  in  same  regj 
on  plate  with  uniform  wall  temperature  0  -  0Q  .  Abscissa  is  value  < 
0/0  at  x;  r  is  an  exponent  in  equation  (17 )  (relation  by  wich  wai: 
temperature  variation  Is  prescribed) j  Prandtl  number  Is  parameter  on 
curves . 
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(b)  r  -  1/2. 

Figure  2.  -  Continued.  Ratio  of  over-all  heat  transfer  from 
x  =  0  to  x  on  plate  with  variable  wall  temperature  to  over¬ 
all  heat  transfer  in  same  region  on  plate  with  uniform  wall 
temperature  9  ■=  @Q  .  Abscissa  is  value  of  0/9o  at  r  ls 
an  exponent  is  equation  (17)  (relation  by  which  wall-tempera¬ 
ture  variation  is  prescribed);  Prandtl  number  is  parameter  on 
curves . 
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Figure  2.  -  Continued.  Ratio  of  over-all  heat  trans- 
fer  from  x  =  0  to  x  on  plate  with  variable  wall 
temperature  to  over-all  heat  transfer  in  same  re¬ 
gion  on  plate  with  uniform  wall  temperature  0  =  6Q 
Abscissa  is  value  of  0/0o  x»  I’  ^-s  311  exponent 

in  equation  (17) (relation  by  which  wall -temperature 
variation  is  prescribed);  Prandtl  number  is  param¬ 
eter  on  curves. 
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(e)  r  =  3. 


Figure  2.  -  Concluded.  Ratio  of  over-all  heat  trans¬ 
fer  from  x  =  0  to  x  on  plate  with  variable  wall' 
temperature  to  over-all  heat  transfer  in  same  re¬ 
gion  on  plate  with  uniform  wall  temperature  (0  =  e  ) 
Abscissa  is  value,  of  0/0o  at  x;  r  is  an  exponent 
in  equation  (17) (relation  by  which  wall -temperature 
variation  is  prescribed);  Prandtl  number  is  param¬ 
eter  on  curves. 
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(a)  r  »  0. 


Figure  3.  -  Ratio  of  local  heat-transfer  coefficient  at  x  for 
plate  with  variable  wall  temperature  to  local  coefficient  at 
same  location  on  plate  with  uniform  rail  temperature  6  c  6 
Abscissa  is  value  of  6/9  0  at  x;  r  is  an  exponent  in  equa¬ 
tion  (17) (relation  by  which  wall-temperature  variation  is  pre¬ 
scribed)  ;  Prandtl  number  is  parameter  on  curves . 
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(b)  r  -  1/2. 

Figure  3.  -  Continued.  Ratio  of  local  heat-transfer  coefficient 
at  x  for  plate  with  variable  wall  temperature  to  local  coeffl 
cient  at  same  location  on  plate  with  uniform  wall  temperature 
«  ^  - - - a/o  of  v  r*  is  an  exponent 


«  a  «n  .  Abscissa  is  value  of  w/vo  *“* 

in  equation  (17) (relation  by  which  wall-temperature  variation 
is  prescribed);  Prandtl  number  is  parameter  on  curves. 
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Figure  3.  -  Continued.  Ratio  of  local  heat-transfer  coefficient 
at  x  for  plate  with  variable  wall  temperature  to  local  coeffi¬ 
cient  at  same  location  on  plate  with  uniform  wall  temperature 
0  =■  60  .  Abscissa  is  value  of  0/0  at  x;  r  is  an  exponent 
in  equation  ( 17) (relation  by  which  wall-temperature  variation 
is  prescribed) ;  Prandtl'  number  is  parameter  on  curves . 


mud 
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Figure  3.  -  Concluded.  Ratio  of  local  heat-transfer  coefficient 
at  x  for  plate  with  variable  wall  temperature  to  local  coeffi¬ 
cient  at  same  location  on  plate  with  uniform  wall  temperature 
8  -  0  .  Abscissa  is  value  of  0/0o  at  x;  r  is  an  exponent 

in  equation  (17 ) (relation  by  which  wall-temperature  variation 
is  prescribed);  Prandtl  number  is  parameter  on  curves. 
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